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Abstract 

We present a new probabilistic algorithm to find a finite set of points intersecting 
the closure of each connected component of the realization of every sign condition over 
a family of real polynomials defining regular hypcrsurfaces that intersect transversally. 
This enables us to show a probabilistic procedure to list all feasible sign conditions 
over the polynomials. In addition, wc extend these results to the case of closed sign 
conditions over an arbitrary family of real multivariate polynomials. The complexity 
bounds for these procedures improve the known ones. 

1 Introduction 

Given polynomials /i, • • • , f m £ M.[x\, . . . , x n ], a sign condition cr G {<, = , >} m , or a closed 
sign condition a G {<, =, >} m , is said to be feasible if the system /i(x)<ti0, • • • , / m (i)a m 
has a solution in M n , and the set of its solutions is called the realization of a. One of 
the basic problems in computational semialgebraic geometry is to decide whether a sign 
condition is feasible. This problem is a particular case of quantifier elimination and, on 
the other hand, many elimination algorithms use subroutines determining all the feasible 
sign conditions for a family of polynomials. 

The first elimination algorithms over the reals are due to Tarski [¥T] and Seidenberg 
[39] , but their complexities are not elementary recursive. Collins [16] was the first to obtain 
a doubly exponential complexity. In [21], Grigor'ev and Vorobjov present an algorithm 
with single exponential complexity to decide the consistency of a system of equalities 
and inequalities by studying the critical points of a function in order to obtain a finite 
set of points intersecting each connected component of the solution set. This same idea 
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was used to obtain more efficient quantifier elimination procedures (see [25], [32j and 
[7])- A standard technique is to take sums of squares and introduce infinitesimals to 
reduce the problem to the study of a smooth and compact real hypersurface. The specific 
problem of consistency for equalities over R was also treated through the critical point 
method afterwards. In [33] the non-emptiness of a real variety defined by a single equation 
is studied, reducing the introduction of infinitesimals and, in [2], an algorithm with no 
infinitesimals is given to deal with arbitrary positive dimensional systems. 

Several probabilistic procedures lead to successive complexity improvements. Using 
classical polar varieties, in [3] and [3], the case of a smooth compact variety given by a 
regular sequence is tackled within a complexity depending polynomially on an intrinsic 
degree of the systems involved and the input length. To achieve this complexity, straight- 
line program encoding of polynomials and an efficient procedure to solve polynomial equa- 
tion systems over the complex numbers ([IE]) are use d- The compactness assumption is 
dropped in [5] and [6], by introducing generalized polar varieties. The non-compact case 
is also considered in [36j for a smooth equidimensional variety defined by a radical ideal 
by studying projections over polar varieties, and an extension to the non-equidimensional 
situation is presented in [37]. Finally, [35] deals with sets of the type {/ > 0} through the 
computation of generalized critical points. 

In this paper, we consider the problem of determining all feasible sign conditions (or 
closed sign conditions) over a given finite family of multivariate polynomials. We first 
present a probabilistic algorithm that, under certain regularity conditions (see Hypothesis 
[9]), obtains a finite set of points intersecting the closure of each connected component of 
the realization of every sign condition over the given polynomials. For families of arbi- 
trary polynomials, we show a probabilistic algorithm that computes a finite set of points 
intersecting each connected component of the realization of every closed sign condition. 
The input and intermediate computations in our algorithms are encoded by straight-line 
programs (see Section [2.2p . The output is described by means of geometric resolutions, 
that is to say, by univariate rational parametrizations of 0-dimensional varieties. In both 
situations, the output of the algorithm enables us to determine all the feasible closed sign 
conditions over the polynomials by evaluating their signs at the computed points, which 
is done by using the techniques in [T5]. Moreover, in the first case, we can determine all 
feasible sign conditions (for a treatment of the problem in full generality, see for instance 
0). 

A sketch of our main algorithms is the following. Given . . . , f m £ R[xi, . . . , x^], a 
generic change of variables prevents asymptotic behavior with respect to the projection 
to the first coordinate x\ for each connected component C C R n of every feasible (closed) 
sign condition over fi, . . . , f m : either C projects onto R or its projection is a proper 
interval whose endpoints have a nonempty finite fiber in C. In the latter case, points in 
C are obtained as extremal points of x\. These extremal points are solutions of particular 
systems of polynomial equations which are dealt with by deformation techniques that 
enable the computation of geometric resolutions of finite sets including them. To find 
points in the components projecting onto R, the set is intersected with {x\ = pi} for a 
particular value pi, and the algorithm continues recursively. 

The following theorem states our main results (see Theorems \W\ and [26]) : 
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Theorem Let K be an effective subfield of K. There are probabilistic algorithms to 
perform the following tasks: 

• Given polynomials fi, ■ ■ ■ , f m £ . . . , x n ] satisfying Hypothesis [3 mi/i degrees 
bounded by d > 2 and encoded by a straight-line program of length L, obtain a finite 
set of points intersecting the closure of each connected component of the realiza- 
tion of every sign condition over fi, . . . , f m within 0( ^™^{ m > n } (("Zi)^ n ) 2 -^) 
operations in K up to poly-logarithmic factors. 

• Given arbitrary polynomials fi, ■ ■ ■ , f m € K[x\, . . . , x n ], with degrees bounded by an 
even integer d and encoded by a straight-line program of length L, obtain a finite set 
of points intersecting each connected component of the realization of every closed sign 
condition over fi, ■ ■ ■ , f m within 0( Y^^ m ' n ^ 2 s (™) (("~J)d n ) 2 (L + d)) operations 
in IK up to poly-logarithmic factors. 

The factor (™Zi)^ n i n the complexity estimates is an upper bound for the bihomo- 
geneous Bezout numbers arising from the Lagrange characterization of critical points of 
projections (cf. [37J). In fact, one of the new tools to achieve the stated complexity or- 
der, which improves the previous ones depending on the same parameters, is the use of 
algorithmic deformation techniques specially designed for bihomogeneous systems (for a 
similar approach, see [23] ). Up to now, the polynomial systems used to characterize crit- 
ical points were handled with general algorithms solving polynomial equations over the 
complex numbers (see, for instance, [I], [33], [20] and [30]). As taking sums of squares of 
polynomials and introducing infinitesimals lead to an artificial growth of the parameters 
involved in the complexity estimates, another important feature of our techniques is that 
we work directly with the input equations instead of using these constructions. 

This work can be seen as an extension of [36] and of [6| in the sense that we deal not 
only with equations but also with inequalities. In particular, the algorithm in [36] . which 
only works for the case of smooth equidimensional varieties defined by a radical ideal, 
considers a family of equation systems equivalent to the ones introduced in the recursive 
stages of our algorithm. However, those systems involve a large number of polynomials 
and do not have any evident structure. 

We also prove that our deformation based approach can be applied to deal with sign 
conditions over bivariate systems without any assumption on the polynomials. We expect 
this can be extended to general multivariate polynomials. This is the subject of our 
current research. Finally, we adapt our techniques to the case of an arbitrary multivariate 
polynomial. 

This paper is organized as follows: In Section [2] we introduce some basic notions and 
notation that will be used throughout the paper. Section [3] is devoted to presenting the 
basic ingredients to be used in the design of our algorithms. In Section IU we present 
our main algorithms to determine all feasible sign conditions over polynomial families 
satisfying regularity assumptions. In Section [5] we consider the same problem for closed 
sign conditions over arbitrary multivariate polynomials. The last section contains our 
results on sign conditions over bivariate polynomial families and over a single multivariate 
polynomial. 



3 



2 Preliminaries 



2.1 Notation 

Throughout this paper Q, K and C denote the fields of rational, real and complex numbers 
respectively, N denotes the set of positive integers and No := NU {0}. If A; is a field, k will 
denote an algebraic closure of k. 

For n G N and an algebraically closed field k, we denote by A£ and (or simply by 
A n or P n if the base field is clear from the context) the n-dimensional affine space and 
projective space over k respectively, equipped with their Zariski topologies. For a subset 
X of one of these spaces, we will denote by X its closure. We adopt the usual notions of 
dimension and degree of an algebraic variety V (see for instance |40j and [22j). 

We will denote projections on certain coordinates x by tt x . For short, a projection on 
the kth coordinate will also be denoted by TTk- 

For any non-empty set A C W 1 , A will denote its closure with respect to the usual 
Euclidean topology. We define Z in f(A) = {(x\, . . . ,x n ) G A | x\ = inf tti(A)} if -k\(A) is 
bounded from below, and Z m {(A) = otherwise. Similarly, Z sup (A) = {(cci, . . . ,x n ) E A \ 
x\ = sup7Ti(yl)} if tti(A) is bounded from above, and Z snp (A) = otherwise. Finally, we 
denote Z(A) = Z ini (A) U Z sup (A). 

Throughout this paper, log will denote logarithm to the base 2. 

2.2 Algorithms and complexity 

The algorithms we consider in this paper are described by arithmetic networks over an 
effective base field KcK (see [33] ) . The notion of complexity of an algorithm we consider 
is the number of operations and comparisons over K. 

The objects we deal with are polynomials with coefficients in EC. Throughout our 
algorithms we represent each polynomial either as the array of all its coefficients in a 
pre- fixed order of its monomials (dense form) or by a straight-line program. Roughly 
speaking, a straight-line program (or sip, for short) over K encoding a list of polynomials 
in K[sei, . . . , x n ] is a program without branches (an arithmetic circuit) which enables us to 
evaluate these polynomials at any given point in W 1 . The number of instructions in the 
program is called the length of the sip (for a precise definition we refer to [141 Definition 
4.2]; see also [26]). 

We will do operations with polynomials encoded in both these ways. To estimate the 
complexities we will use the following results: Operations between univariate polynomials 
with coefficients in a field EC of degree bounded by d in dense form can be done using 
0(d\og(d) log log(d)) operations in EC (see [441 Chapter 8]). From an sip of length L en- 
coding a polynomial / G K[x±, . . . , x n ], we can compute an sip with length O(L) encoding 
/ and all its first order partial derivatives (see [10j). 

2.3 Geometric resolutions 

A way of representing zero-dimensional affine varieties which is widely used in computer 
algebra nowadays is a geometric resolution. This notion was first introduced in the works 
of Kronecker and Konig in the last years of the XlXth century ([29] and [28]) and appears 
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in the literature under different names (rational univariate representation, shape lemma, 
etc.). For a detailed historical account on its application in the algorithmic framework, 
we refer the reader to [20]. The precise definition we are going to use is the following: 

Let fcbea field of characteristic and V = . . . , £ (D) } C A£ be a zero-dimensional 
variety defined by polynomials in k[x\, . . . , x n ]. Given a separating linear form I = u\X\ + 
• • - + u n x n G k[x\, . . . ,x n ] for V (that is, a linear form t such that / l{$>)) if i^j), 

the following polynomials completely characterize the variety V: 

• the minimal polynomial q := Tli<i<D(^ — ^(C®)) £ M^l °f ^ over the variety V 
(where U is a new variable), 

• a polynomial q G &[£/] with deg(q) < D and relatively prime to q, 

• polynomials w\, . . . , w n G k[U] with deg(u»,) < -D for every 1 < j < n satisfying 

V={Q( V ),...,'^(r l ))eT\r l ek, q{ n ) = 0}. 

The family of univariate polynomials q,q,w\, . . . ,w n £ k[U] is called a geometric resolution 
of V (associated with the linear form £). 

We point out that the polynomial q appearing in the above definition is invertible in 
k[U]/(q(U)). Setting Vk(U) := q~ 1 (U)utk(U) mod (q{U)) for every 1 < k < n, we are lead 
to the standard notion of geometric resolution: a family of n + 1 polynomials q,v±, . . . ,v n 
in k[U] satisfying V = { (^i(^), . . . , v n {ri)) G k | r/ G A;, q(r/) = 0}. We will use both 
definitions alternatively, since the complexity of passing from one representation to the 
other does not modify the overall complexity of our algorithms. Which notion is used in 
each case will be clear from the number of polynomials. 

3 General approach 

3.1 Avoiding asymptotic situations 

For any non-empty set A C W 1 we define Z- m {(A, k) = {(x±, . . . , x n ) G A \ x^ = inf ^^(A)} 
if A is bounded from below and Z[ n f(A, k) = otherwise. Similarly, Z snp (A, k) = 
{(x\, . . . , x n ) G A | Xk = sup7Tfc(^4)} whenever A is bounded from above and Z sup (A, k) = 
otherwise. Finally, Z(A, k) = Z^(A, k)L)Z sup (A, k). In particular, when k = 1, Z[ n f(A, 1), 
Z sup (A, 1) and Z(A, 1) will be denoted by Z- m t (A), Z sup (A) and Z(A) respectively as has 
already been stated in Section 12.11 

The precise conditions achieved by a generic linear change of variables are stated in 
the following: 

Proposition 1 Let fi,...,f m be n-variate polynomials with real coefficients. After a 
generic linear change of variables over Q, for every semialgebraic set V defined in M. n by 
a Boolean formula on the polynomials /i,...,/ m involving equalities and inequalities to 
zero and every p = (pi, . . . ,p n ) G M n , if 1 < k < n and C is a connected component of 
V n {xi = pi, . . . , Xk_i = Pk~i}, then Z(C, k) is a finite set (possibly empty). Moreover, if 
T^kiC) is bounded from below, then Z- m i(C,k) is not empty, and, if vrfc(C) is bounded from 
above, then Z sup (C,k) is not empty. 
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To prove Proposition Q3 we will use the following auxiliary lemma: 

Lemma 2 Let {fij}i<i< n ,i<j<ii C M[xi, . . . , x n ] be a family of nonzero polynomials satis- 
fying simultaneously: 

a) for 1 < i < n, {fij}i<j<ii is contained in R[xi, . . . ,Xj], it is closed under derivation 
with respect to the variable x i} and every polynomial in it is quasi-monic (that is, 
monic up to a constant) with respect to Xi, 

b) for 1 < i < n, {f(i-i)j}i<j<l i _ 1 slices {fij}i<j<k in the sense of LH Definition 
2.34]. 

Let p = (pi, . . . ,p n ) G W 1 , 1 < % < n, and P C i' be a semialgebraic set defined by a 
Boolean formula on the polynomials fij, 1 < j <h, involving equalities and inequalities to 
zero. For 1 < k < i, let C be a connected component ofVD {x\ = pi, . . . = Pk-i}- 

Then the set Z(C,k) is finite (possibly empty). Moreover, ifir^C) is bounded from below, 
then Z m {{C,k) / 0, and, if irk{C) is bounded from above, then Z snp (C,k) ^ 0. 

Proof: As for every 1 < k < n the family {fijipi,. ■ ■ ,Pk-i,x k , ■ ■ ■ , Xn)}k<i<n,i<j<ii c 
M[xfc, . . . , x n ] satisfies the hypotheses, it is enough to prove the lemma for k = 1. 
For i = 1, the result is clear. 

Suppose the statement is true for i — 1. Let 7r : W — > W" 1 be the projection on 
the first i — \ coordinates. Following the notation in |13^ Chapter 2], let A\, . . . ,A# be 
the semialgebraic sets giving the slicing of with respect to fn, . . . , f%. given by the 
polynomials . . . , f\ i -i)i i _ 1 and, for 1 < s < £, let £ Sjl < ■ ■ ■ < £ S)Cls : A s R be the 

continuous semialgebraic functions that slice A s x M. Let A s ^i, . . . , A StUg be the connected 
components of A s . 

Note that C is a finite union of some sets of the partitions of the sets A S}U x M given by 
£ Sj i, . . . ,£, s ,a a and 7r(C) is a finite union of sets A S;U . If 7r(C) = UhA Sh ^ Uh , then Z(n(C)) C 
UhZ(A ShtUh ). Since each A ShtUh is a connected component of A Sh , which can be described 
by a Boolean formula involving equalities and inequalities to zero of fu-i)i, ■ ■ ■ , /(i— ; 
by inductive hypothesis, each Z(A ShMh ) is finite and, therefore, Z(ir(C)) is finite too. 

Let w e Z(C). As Z(C) C C, vr(u;) E vr(C) C vr(C). Moreover, as vri(C) = ^i(^(C)), 
7r(it;) G Z(k(C)). On the other hand, as w £ Z(C), at least one of the quasi-monic 
polynomials fa, - ■ ■ , fn i vanishes at w (otherwise, all the fij, j = 1, . . . , k, would have a 
constant sign in a neighborhood of w contained in C). Therefore, Z{C) is a finite set 
since tt{w) lies in the finite set Z(tt(C)) and Wi is a zero of one of the nonzero polynomials 
fa(7r(w),Xi), . . . ,f ih (7r(w),Xi). 

Suppose now that vri(C) is an interval bounded, for example, from below. Then, by 
the inductive assumption, there exists z = (z±,... , Zi-\) G Zi n f(7r(C)) C vr(C). Assume 
further that A\^ C vr(C), z G Z in { (Ai t i) and 7 : [0, 1] — > Ai t i is a continuous semialgebraic 
curve such that 7((0, 1]) C A\^ and 7(0) = z (see [T31 Theorem 2.5.5]). Let x = 7(1). 
Since x G A\^ C vr(C), there exists i/El such that G C. Using [T3l Lema 2.5.6], 

each £i t a can be extended continuously to A\. Let us denote by ^i )0 also this extension. 
Depending on the position of y with respect to the values £1,1(2;) < • • • < £i, ai (x), it is 
easy in any case to define a continuous semialgebraic function h : [0, 1] — * R such that the 
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continuous function 7 : [0, 1] — ► W denned as j(t) = (7(i), h(t)) satisfies 7((0, 1]) C C (note 
that the signs of the polynomials fn, . . . , fu t are constant over 7((0, 1])) and, therefore, 
(z, h(0)) = 7(0) G C. Moreover, as z x = inf 7ri(7r(C)) = inf vri(C), (z, /i(0)) G Z inf (C). □ 

Now, we can prove Proposition [1} 

Proof: By Lemma [21 it suffices to show that there exists a Zarisky open set U C Gl(n, C) 
such that, for every Vq G Q nxn r\U there exists a family of polynomials {fij}i<i<n,i<j<li c 
satisfying the hypotheses of the lemma, and such that, for 1 < j < m, f n j( x ) — 
/j(Vox) with m < / n . Let V be a matrix whose entries are new variables v rs , 1 < r, s < n 
and consider {Fij}i<i<n,i<j<k C defined in the following way: 

• Take l' n = m and, for 1 < j < l' n , let F n j(V,x) = fjiVx). Then, for 1 < jo < l' n , if 
deg^ F n j = d n j , add the first d n j — 1 derivatives of F n j with respect to x n to the 
list to obtain {F n j}i<j<i n - 

• Prom{F (io+1)j } 1 < :7 <i. o+1 C R[v,xi, ... ,x io+ i], first, form {F^}^^ cK^ii,...,^] 
by taking all possible resultants and subresultants with respect to the variable Xj +i 
between pairs of polynomials, not taking into account the ones that are identically 
zero. Then, for 1 < jo < l' io , if deg x Fi j = di j , add the first di j — 1 derivatives 

of Fi Q j Q with respect to the variable Xi to obtain the family {-Fi j}i<j<z io ■ 

Let 1 < i < n and 1 < j < h- Let dij := deg^ Fij and let qij G M.[v] be the coefficient 
of the monomial x^ J in Fij G We will show that q^ ^ 0. To do so, we also 

prove the following: for A G Q JXJ , denote by A and A the matrices A = f ^ ^ ~\ G 

Q(i+i)x(i+i) an d i = ( J G Q nxn ; then, for every I < j < l[ and every 

5 G Q ixi , Fij(V,Bx) = Fij(VB,x), and, for every Z J + 1 < j < ^ and B G QC*- 1 )*^- 1 ), 
FyCV.Ss) =F ii (y J B,x). 

If i = n, it is clear that q n j 7^ for every 1 < j < l' n because there is a change of 
variables Vo G Q nxn which makes the polynomials fj, for 1 < j < m, quasi-monic of 
degree d n j with respect to the variable x n . Moreover, for 1 < j < l' n and for B G <Q) nxri , 
F nj {V,Bx) = fj(VBx) = F nj {VB,x). For l' n + 1 < j < l n and B G qC™- 1 )*^" 1 ), 
suppose F nj - is the derivative with respect to x n of F n j/(V,x). Then, by differentiating 
F n j'(V, Bx) = F n ji(VB,x) with respect to x n , we obtain the desired identity. 

Let i < n and B G Q ixi . Note that, if 1 < j < and 

d (i+l)j 



F(i+i)j(V,x) = ^2 g( i+1 )^(y,xi,... 



h=0 



as = F {i+1)j (VB,x), then g (i+1 ) i)h (V, Bx) = g( i+1 ) i)/l (V r -B,a;) for 1 < h < 

da + xy. Let 1 < j < and suppose F^ is a resultant or a subresultant between i^j+i)! 
and F {i+1)2 . As deg x . +1 F {i+1)1 = d {i+1)1 and deg x . +1 F (i+1 ) 2 = d(»+i)2, corresponds 
to the evaluation of a certain polynomial i? in qu+uih and Q r (i+i)2,/i- Then Fij{V,Bx) = 



7 



R(l{i+i)i,l(V,Bx),.. . ,1(i+i)2,d (i+1)2 (V,Bx)) = R(q( i+ i)i,i{VB,x), . . . , g(i+i)2,d Ci+1)2 (VB, x)) = 

Fij(VB,x). When + 1 < j < Zj, we may proceed in the same way as we do for 
i = n. Now, it suffices to prove that qij ^ for 1 < j < l^. Let Bo £ Q IXI such 
that deg Xi Fij(V,B x) = d^. As Fij(V,B x) = Fij(VB Q ,x), then qij(VB ) / and, 
therefore qij ^ 0. 

Define U = {Vq £ C nxn \ qij{V ) / for 1 < i < n, 1 < j < k}. By Proposition 
4.34 and Theorem 5.14], for every V £ Q nxn nZ^, the set {/jj(x)}i<j< ni i<j<; i defined by 
fij(x) = Fij(Vo,x) satisfies both conditions in Lemma [2j □ 

The following proposition is a major tool for our algorithms (cf. [361 Theorem 2]). 

Proposition 3 Let /i,...,/ m be n-variate polynomials with real coefficients. After a 
generic change of variables, for every semialgebraic set V defined in W 1 by a Boolean for- 
mula on fi, . . . , f m involving equalities and inequalities to zero, and every p = (pi, ■ ■ ■ ,p n ) €= 
M. n , if for 1 < k < n, V(k,p) is the set of all the connected components of V n \x\ = 
pi, . . . ,x k -i = Pk-i}, then 

n 

Mu(U |J z(c,k)) 

k=l C&V(k,p) 

is finite and intersects the closure of each connected component ofV. 

Proof: Proposition Q] ensures that this set is finite. Let C\ £ V(l,p). If 7Ti(Ci) is 
bounded from above or below, again Proposition Q] states that Z(C\,1) is a finite non- 
empty set and is included in C\. Otherwise, 7Ti(Ci) = E and C\ n {x\ = Pi} ^ 0. Let 
C2 £ V(2,p) be a connected component of C\ n \x\ = p{\. If ^(C^) is bounded from 
above or below, Z(C2,2) / and is included in C2 C Ci. Otherwise, ^2(62) = R and 
Ci fi {xi = pi, X2 = P2} 7^ 0- Following this procedure, we obtain that either there exists 
Ct £ V(k,p) such that Z(Ck,k) 7^ and is included in Ck C Ci for some 1 < k < n or 
p£Ci. □ 

The proof above leads to the recursive structure of our algorithm: For 2 < k < n, we 
may think the fc-th variable as the first one for the polynomials fj(pi, ■ ■ ■ ,pk-i,Xk, ■ ■ ■ ,x n ) 
for 1 < j < m. Therefore, it is enough to consider the problem of finding extremal points 
for the projection over the first coordinate of the closures of the connected components of 
a semialgebraic set. 



3.2 Equations defining extremal points 



Let /1, . . . , f m £ M[xx, . . . , x n ] and let 5 := {i\, . . . , i s } C {1, . . . , m}. If 1 < s < n — 1, 
the implicit function theorem implies that a point z with maximum or minimum first 

0} satisfies 



coordinate in a connected component of {f - tl = ■ ■ ■ = fi, 
fh( z ) =■■■ = fi s (z) = 0, rank 



V !&(*) 



dx n \ Z > \ 



< s. 



(1) 
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This condition can be rewritten as 



f h (z) =^. = f ls ( z ) = 0, 



0). 



(2) 



for fj,x, . . . , fj, s G M not simultaneously zero, where V/j -(z) denotes the vector obtained by 
removing the first coordinate from the gradient \/fij(z). 
When s = 1, for z G M n , conditions (pQ) are equivalent to 




When s > ?i, we will simply consider the conditions 



fM = o. 



(4) 



For every 2 < s < n — 1, as system ([2]) is homogeneous in the variables /xi, . . . , fi s , we 
consider the variety Ws C AJi x P^" 1 defined as the zero set of this system. If s = 1 or 
s > n, let Ws be the variety defined by systems ([3]) and (jlj) respectively. 

The following result is an adaptation to our context of the Karush-Kuhn- Tucker con- 
ditions (see [31J) from non-linear optimization which generalize the Lagrange multipliers 
theorem in order to consider equality and inequality constraints. 

Proposition 4 Let /i, . . . , f m € R[xi, . . . ,x n ] and a = (a x , . . . ,a m ) G {<, <,=, >, >} m . 
Set E a = {i | (Tj = " = "}. Then, for every connected component C of the set {x G M n | 
/i(x)<TiO, . . . , /m^ffmO}, we have 



Proof: Without loss of generality, assume E a = {1,. . . ,1} (or E a = and I = 0) and 
7Ti(C) is bounded from below. Let z = (zi, . . . , z n ) € Z m f(C) and So = {« € {1, . . . , m} | 
/i(^) = 0}. Note that C So and, even when I = 0, So 7^ 0; then, we may assume that 
So = {1, . . . , t} with max{l, 1} < t < m. We will show that z G tTx(Ws )- 
If t > n, we have z G ir x (Ws ) by the definition of this set. 

Assume now that t < n — 1. If z ^ tt x (W5 ), the set {Vfi(z),i G So} is linearly 
independent. Let / : M n — > K* be the map / = (/1, . . . , ft). We may assume that the minor 
corresponding to the variables n — t + 1, . . . , n in the Jacobian matrix Df(z) is not zero. 
Applying the inverse function theorem to h(x) = (x±— z\, . . . , x n -t — z n _ t , fi(x), ■ ■ ■ , ft(x)), 
there exist an open neighborhood U of z, e G M>o and a map g : (— e, e) n — > U inverse to 
h : U — ► (— e, e) n . Moreover, we may assume that /t+i, . . . , f m have constant signs on [/. 

Let w £ CflP and let y = h(w). Let a G {<,=,>} m be such that fi(w)ai0 for 
1 < i < m. Then, the conditions y\ = w\ — Z\ > 0,y n -t+iai0, . . . ,y n at0 hold. Since 
w G C, for 1 < i < m, oi G {<, =} if Oi = " < ", a« G {>, =} if <Jj = " > " and 0{ = Oi 
otherwise. Hence, every point satisfying a also satisfies a. 

Let 7 : [-e/2,yt] -> (-e,e) n be defined as 7(1*) = (n, y 2 , . . . ,y n ). For u G [-e/2,yi] 
and 1 < i < t, o 7(n)) = y n „ i+ j(7i0. Taking into account that, for t + 1 < i < m, fi 



Z(C) c |J 



B CT CSc{l,...,m}, 
S^0 
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has constant sign over U and the image of goj lies in U, we also have that fi{g°y{u))ai{) 
for t + 1 < i < m. Therefore, the image of g o 7 is contained in the realization of a 
and, since it is a connected curve with a point g o j(y\) = w in the connected component 
C, we conclude that it is contained in C. Now, the first coordinate of g o y(—e/2) is 
— e/2 + z\ < z\, contradicting the fact that z\ = inf iti(C). □ 

3.3 Deformation techniques for bihomogeneous systems 

In this subsection we present briefly a symbolic deformation introduced in [19] . |18j . |24j . 
|20j and [38], adapted to the bihomogeneous setting following [23]. 

3.3.1 The deformation 

Given polynomials h\ (a;), ... , h s (x), h s+ i(x, fi), . . . , h r (x, fj) E K[x±, . . . , x n , m, . . . , fj, s ] we 
consider the associated equation system: 

hx(x) = 0, . . . , h s {x) = Q,h s+1 (x,(j,) = 0, . . .,h r (x,fi) = 0, (5) 

where 2 < s < n — 1 and r = s + n — 1 such that, for 1 < i < r, deg x (/ij) < di < d and, 
for s + 1 < i < r, hi is homogeneous of degree 1 in the variables [i. 

Let W C A n x P s_1 be the variety this system defines. By the multihomogeneous 
Bezout theorem ([3UJ Ch. 4, Sec. 2.1]), the degree of W is bounded by 

i=l Bc{s+l,...,r},#fi=n-sj6S v 7 

Let 51 (x), . . . ,g s (x),g s+ i(x,n), . . .,g r (x,n) £ K[xt, . . . ,x n , m, . . .,/i s ] be polynomials 
with deg x (gi) = di for 1 < i < r and homogeneous of degree 1 in the variables \x for 
s + 1 < i < r, such that: 

(H) gi , . . . , g r define a 0-dimensional variety in A n x {fi s ^ 0} C A n x P s_1 with D points 
si, . . . , so satisfying 7r x (sj) 7^ n x (sj) for i 7^ j, and the Jacobian determinant of the 
polynomials obtained from gi, . . . , g r by dehomogenizing them with (j, s = 1 does not 
vanish at any of these points. 

We will specify polynomial systems meeting these conditions in Definitions H2] and [201 
below. 

Let t be a new variable. For every 1 < i < r, let 

Fi := (1 - t)hi + t 9i . (7) 

Consider the variety V C A 1 x A n x F s_1 defined by Fi, . . . , F r , and write 

V = V {0) U V {1) U V, (8) 

where is the union of the irreducible components of V contained in {t = 0}, is 
the union of its irreducible components contained in {t = to} for some to G C \ {0}, and 
V is the union of the remaining irreducible components of V. 
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Lemma 5 With our previous assumptions and notation, tt x u{V Pi {t = 0}) is a finite 
subset of W containing all its isolated points. 

Proof: Let V\ be an irreducible component of V and V\ be its Zariski closure in A 1 x P n x 
P s_1 . The projection of V\ to A 1 is onto and so, V\ D {t = 1} 7^ 0. But our assumption 
on gi,...,g r implies that Vi n {t = 1} = V\ n {t = 1}; then, Vi n {t = 1} 7^ and it is 
0-dimensional. It follows that dim(Vi) = 1. Therefore, V is a 1-equidimensional variety 
and thus n X)IJ ,(V n {t = 0}) is a finite set. 

In order to prove the second part of the statement, note that W = ir x ^(V n {t = 0}) = 
^x,fi(V^) U TTx^iV n {t = 0}). Now, an isolated point of W cannot belong to 7r X)M (V(°)), 
since the dimension of each of its irreducible components is at least 1; hence, it lies in 

7Tx lM (^n{t = 0}). □ 

The same deformation can be applied to a system h\(x), . . . , h n (x) G K[x\, . . . ,x n ] 
with gi (x), . . . , g n {x) G . . . , x n ] such that deg(^) = deg(/tj) for every 1 < i < n and 

having Y\a=\ deg(gi) common zeros in A n . 

3.3.2 A geometric resolution 

Lemma 6 The variety defined in A^j x P^r^y by Fi, . . . , F r as in (Q) is 0-dimensional 

and has D points S\, . . . , Sd in {fj, s 7^ 0} such that 7r x (Si) 7^ ir x (Sj) for i 7^ j. Moreover, 
these points can be considered as elements in — l]] r . 

Proof: The multihomogeneous Bezout Theorem (see, for instance, |4CH Chapter 4, Section 
2.1]) states that the degree of the variety is bounded by D. If Sj, 1 < i < D, are 
the common zeros of g%, ... , g s ,g s +\, . . . , g r , the Jacobian of F\, . . . ,F r with respect to 
zi, . . . , x n ,{j,±, . . . , Hs-i at t = 1 and (x, /x) = Sj is nonzero. The result follows applying 
the Newton-Hensel lifting (see for example [241 Lema 3]). □ 

Consider now new variables y±, . . . , y n and define £(x, fx, y) = i{x, y) = Yuj=i Vj x j- F° r 
ai, . . . , a n e C, let £ a (x, y) = l a (x) = YTj=i a j x j- Let 

P(t,U,y) = n (U-t{Si, V )) = E ° =0 l h { ^ y)Uh = ^^f 1 e K(t)p7,y], (9) 

with P(t, J7, y) G f7, y] with no factors in K[t] \ K. 

In order to compute P we will approximate its roots. The required precision is obtained 
from the following upper bound for the degree of its coefficients. A similar result in the 
general sparse setting appears in [27} Lemma 2.3], but to avoid a possibly cumbersome 
translation of our setting into sparse systems, we give an alternative statement and its 
proof here. 

Lemma 7 Using the notation in deg t P(t, U, y) < nD. 
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Proof: Let : V x A n — ► A n+2 be the morphism defined by $>(t,x, n,y) = (t,£(x,y),y). 
It is easy to see that P(t, U, y) is a square-free polynomial defining Im<3?. 

For = (f3 , AO E C n +\ let P^(t) = P(t, A, A, • • • , AO- For a generic A 

Pg(t) is squarefree and satisfies deg t P(t, U, y) = deg t Pp(t) = #(lm$ n {£/ = A) y\ = 0i, 
...,y n = At}) - 

Since Im<I> \ Im<I> is contained in an n-dimensional variety, for a generic /3 we have 
Iin¥n {17 = A), yi = A, • • • , Vn = AJ = Im<S>n{U = A, Vi = A, • • • , Vn = AJ- Moreover, 
for a fixed to E C, for a generic (3 we also have Im<I>n {t = to} H {£/ = A> ?/l — A? • • • ? Vn = 
Ai} = 0- Then, for a finite set T C C, for a generic A Im$n{C7 = A)> 2/1 = Aj • • • > JM = Ai} 

= lm$ n ( n toeT {< / to}) n {£/ = A, yi = A, • • • , y n = A*}- 

Let T = {to E C | V has an irreducible component included in {t = to}} and (3 generic 
satisfying all the previous conditions. 

For each (t, (3) E Im$ n (n to£ T{i / *o}) H {17 = A, yi_= A, • • • ,Vn = Ai}> there exists 
a point (t, x,/2) E V such that ^( / g 1| ... i/ a Jl )(^) = A- Then, (t, x, /2) is a solution to the system 

Fi(t, x) = 0, . . . , F s (t, x) = 0, F s+1 (t, x,fi)=0,..., F r (t, x, = 0, %,... i/3n) (x) - A = 

(10) 

with t ^ T. Since P^ is a nonzero polynomial, V has no irreducible components included 
in {£(p 1 .... ) i3 n )(x) = A} an d, therefore, V D {•^( J a 1: ... ) /3„) (a^) = A} is a 0-dimensional variety. 
This implies that the point (t, x,/2) is an isolated solution of the system (fTUj) . Then, in 
order to get an upper bound for deg t P(t, U, y) it suffices to bound the number of isolated 
roots of this system. 

Using the multihomogeneous Bezout theorem in the three groups of variables x, \x 
and t, it follows that this system has at most Yl(E t e x e ) Ylk&E x dk isolated roots, where 
d r+ i := 1 and (E t , E x , E^) runs over all the partitions of {1, . . . , r + 1} into sets E t , E x 
and Efj, of cardinality l,n and s — 1 respectively, with E t C {1, . . . , r} and P M C {s + 
1, . . . ,r}. Each of these partitions (E t , E x , E^) can be obtained uniquely from a subset 
E C {s + 1, . . . , r} of cardinality n — s and an element e E {1, . . . , s} U E by taking 
£ t = { e }, E x = ({1, . . . , s, r + 1} U E) \ {e} and E^ = {s + 1, . . . , r} \ E. Therefore, 

E II dk = J2 II d k < J2 H d k < nD. 

{Et,E x ,Efj) k&E x (E,e) fce({l,...,s,r+l}UE)\{e} (E,e) ke{l,...,s}UE 

□ 

Note that, if tt x (V n {t = 0}) = {zi, . . . , z„} C A n , for every 1 < / < u, U - t(z h y) 
divides P{0,U,y). Then, if P(0,U,y) = Hf =1 qj(U,y) 5 i is the factorization of P(0,U,y) 
in C[U, y], we may suppose that qi = U — £(zi, y) for every 1 < I < v < a. Let Q(U, y) = 

Let a = (ai, . . . ,a n ) E C n such that, for 1 < j < a, qj(U,a) has the same degree 
in U as qj(U,y) and is square- free and, for 1 < ji < j2 < a, Qji(U,a) and qj 2 (U,a) are 
relatively prime polinomials in C[?7]. Note that these conditions are met for a generic 

vector a. Then, P qi^ ,0 \ is square- free and vanishes at £(zi, a) for 1 < I < v. In addition, 
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for 1 < I < v, the kth coordinate of (1 < k < n) is the quotient of 

^(0,l(z t ,a),a) 
Q(£(z/,a),a) 

by 

fg(0,l(z f ,a),a) 
Q(^(z;,a),a) 

We deduce the following: 
Proposition 8 Let P(t,U,y) be as in fSjj. Then, for a generic a G C n , 

l Q(E/,a) ' Q(C/,a) ' Q(C/,a) Q(C/,a) J 

is a geometric resolution of a finite set containing ir x (V n{t = 0}). □ 

4 Regular intersections 

Let fx, ... , f m G R[a;i, . . . , x n ] be polynomials meeting the following condition: 

Hypothesis 9 For every x G C n , for every . . . ,i s } C {1, . . . , no}, if f^ (x) = ■ ■ ■ = 

fi s (x) = 0, then {V/j^x), . . . , V/j s (x)} is linearly independent. 

As in Section 13,21 for S = . . . ,i s } C {1, . . . , m}, consider the solution set Ws of 
the system ([2]), © or Q depending on whether 2 < s < n — 1, s = 1 or s >n respectively. 
Note that, under Hypothesis Ws is the empty set whenever s > n. Moreover: 

Lemma 10 Under Hypothesis 0, after a generic linear change of variables, for every 

5 C {1, . . . , m}, the set Ws is finite. 

Proof: If s = n, Hypothesis [9] implies that Ws is a finite set. 

Now let s < n — 1. Note that n x (Ws) is the set of critical points of the map 
(x±, . . . ,x n ) i— > x\ over the set {x G C n \ fi(x) = for i G S}. By the arguments in 
[421 Section 2.1] based on Sard's theorem and a holomorphic Morse lemma, it follows that 
a generic linear form has a finite number of critical points on this complex variety. There- 
fore, taking any of these generic linear forms as the first coordinate, the set ir x (Ws) turns 
to be finite. Moreover, for every z G tt x (Ws), since {V/j(z), i G S} is linearly independent 
and {Vfi(z), i G S} is linearly dependent, it follows that there is a unique \i G P s_1 such 
that (z,fi) G Ws- □ 

Given /i, . . . , f m G K.[x±, . . . , x n ], we will deal with the deformation (|7|) and the corre- 
sponding varieties ([8]) defined from the systems ([2]), ([3|) and (j4|) for S C {1, . . . , m} with 
1 < #S < n, and an adequate initial system. We will add a subscript 5 in the notation 
V, , V W and V to indicate that the varieties are defined from the polynomial system 
associated with S. 

From Proposition [J Lemma [5] and Lemma [TOl we deduce: 



Si(zi)kY[(lj(£(zi,a),a). 



6iY[qj(£(zi,a),a) / 0. 



13 



Proposition 11 Let fi,---,f m G . . . , x n ] be polynomials satisfying Hypothesis^ 

Let a £ {<,=, >} m and E a = {i | ffj = " = "}. Then, after a generic linear change of 
variables, for each connected component C of {x E M n | /i(x)o"iO, . . . , / m (x)<7 m 0}, 

z(C)c |J 7r*(Vsn{t = o}). 

B (T CSC{l,...,m} 
l<#S<n 

□ 

4.1 Initial systems for deformations 

In what follows we introduce the initial systems for our first algorithmic deformation 
procedure. 

Definition 12 For a given s with 1 < s < n and r := s + n — 1, a type 1 initial system 
is a polynomial system of the form: 

9i{x) = (xi- j) forl<i<s, 

i<j<di 

gi(x,p)=[ Y[ < l>ij( x ))iPi{f J ') fors + l<i<r, 
i<j<di 

where, fors + l<i<r, y (x) = ( £ s+1 < fc < n (i _ a _ 1)d ^_ 1+fc _ 8 gfc) + {i _ s _ 1)d+ }_ 1+n+1 _ s 

(1 < j < di), and ik(fi) = Ei< fc < s i-a-i+fc /^fc- 

For s = 1 and s = n, a type 1 initial system consists of n polynomials of the form 
gi{x) = {xi - j) (l<i<n). 

i<j<d t 

Lemma 13 Property (H) (see Section \3. 3. 1\) holds for any type 1 initial system. 

Proof Let g s +i, ■ ■ ■ ,g r be the polynomials obtained by dehomogeneizing g s +\, . . . ,g r in 
the variables fi with fj, s = 1 and let (xi, . . . ,x n ,p,i, . . . ,fl s ~i) G A™ +s_1 be a solution of 
gi = ■ ■ ■ = g s = gs+i = ■ • • = gr = 0. As, for 1 < i < s, Xi is one of the di zeros of <?j, to 
count the solutions of the original system, the ni<j<s^ solutions (x±, . . . ,x s ) should be 
combined with the solutions . . . , x n , pi, . . . , p s -i) of the system g s+ i = ■ ■ ■ = g r = 0. 

Then, it suffices to show that this new system has ^ec{s+i r} #E=n-s WkeE different 
zeros in A" -1 . 

For s + 1 < i < r, g~i is a product of di linear forms (fiij(x s +i, . . . ,x n ) and the linear 
form tpi(pi, ■ ■ ■ := ipi(pi, ■ ■ ■ j A's— 1)1)- Therefore, each (x s +i, . . . , x n , pi, . . . , Ps—i) 

must be a zero of at least one of these di + 1 linear forms for each i. Suppose it is a zero 
of a linear forms <pij{x) and b > r — s — a linear forms 4>i(fi). Then, . . . , x n , 1) is a 

nonzero vector in the kernel of a Cauchy matrix in Q ax ( n ~ s + 1 ) ) anc l therefore, n — s > a. 
On the other hand, (pi, . . . ,p s -i, 1) is a nonzero vector in the kernel of another Cauchy 
matrix in Q bxs , and therefore s — 1 > b. Asa + 6> r — s = n — 1 = (re — s) + (s — 1) > a + 6, 
then a = n — s and 6 = s — 1. Therefore, for s + 1 < i < r, . . . ,x n ,pi, . . . ,p s ~i) is 

a zero of only one of the di + 1 linear forms defining g~i . 
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The number J2ec{s+i r} #E=n-s TlkeE arises from choosing a subset of n — s 
equations of the form gi = in which a linear form corresponding to the variables x 
vanishes and, then, in each equation, which factor vanishes. In the other s — 1 equations, 
the linear form corresponding to the [i variables must vanish. Any system of this type, 
determines only one point in A 71 " 1 because the Cauchy matrices involved are invertible. 

To see that all the solutions obtained in this way are different, it suffices to show that 
two such solutions satisfy (a;^ 1 , . . . , x„ ) / . . . ,Xn )■ Otherwise, different choices of 
the n — s linear forms would give the same solutions, and this is impossible for any 
choice of n — s + 1 of these linear forms defines an invertible Cauchy matrix. 

Finally, note that the Jacobian matrix of the system gi, . . . , g s ,g s+ i, . . . ,g r is a diagonal 
block matrix. Its upper block is an s x s diagonal block with nonzero elements in its 
diagonal. Without loss of generality, we may suppose that (x s +i, . . . ,x n ) is a common 
zero of the linear forms </>( s +i)i, • • • , <fini and that (fii, . . . , fi s -i) is a common zero of the 
linear forms tp n +ii ■ ■ ■ , Vv- Then, the lower block is the product of a diagonal matrix with 

ds+i d n d„+\ d r 

the elements ^ s+ i(/2) 0( s+ ry(x), . . . ,i/; n (fi) JJ (j) n j(x), ] [ <j)( n +i)j(x), Y\. in 

J=2_ 3=2 j=l _ j=l 

its diagonal by a block diagonal matrix whose blocks are the Cauchy matrices formed by 
the coefficients of the linear forms (t>( s +i)i-> ■ ■ ■ j an d tpn+ij ■ ■ ■ ,ipr respectively, which 
are both invertible. □ 



4.2 Symbolic deformation algorithms 

We describe in this section a probabilistic algorithm which computes a geometric resolution 
as in Proposition O 

Proposition 14 There is a probabilistic algorithm that, taking as input polynomials h\, . . . ,h r 
in K[x±, . . . ,x n ,fii, . . . , fi s ] as in $5\) encoded by an sip of length L, obtains a geometric 
resolution of a finite set containing ir x (V D {t = 0}) for a deformation defined from a type 
1 initial system within complexity 0(n 2 D 2 log(-D) log log(-D) (L + log 2 (-D) log log(-D))) . 

Proof: The procedure of this algorithm is standard. The main difference with previous 
known algorithms solving this task (see, for example, [20] or [23]) is that the Newton 
lifting will be done pointwise. 

First step: Form a type 1 initial system of polynomials of the same degree structure as 
hi, . . . , h r and compute the solutions S\, . . . , sd of this system. The computation of each 
solution amounts to solving two square linear systems of size n — s and s — 1 respectively 
with Cauchy matrices, which can be done within a complexity 0(nlog 2 (n)) by means of 
[121 Ch. 2, Algorithm 4.2]). 

Second step: Construct an sip encoding F±, . . . , F r (see ([7])). Since g±, . . . , g r can be 
encoded by an sip of length 0(dn 2 ), the length of this sip can be taken to be L\ = 
L + 0(dn 2 ). Set F for the list of polynomials F\,...F r dehomogenized with /j, s = 1. 
The algorithm computes, for i = 1, . . . , D, elements Si £ K[t] r such that for 1 < k < r, 

(Si - S l ) k 6 (t — l) 2nD+1 K[[t - 1]]: Let s} 0) = s t be a solution of the initial system 
gi,...,g r . By means of the Newton-Hensel operator we define recursively S^ m+1 ^ = sf"^ — 
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DF-\$ m) )F$ m) ) mod (t - l) 2m+1 K[[t - 1]]. For 1 < k < r and m G N , (^ (m+1) ) fe = 
(§l m) ) k mod (t - l) 2m K[[t - 1]], and (S< m) ) fe is a polynomial in t — 1 of degree less than 
2 m . Operations between polynomials of such degree can be done using 0(2 m mlog(m)) 
operations in K. 

We are now going to estimate the complexity of computing Si '■= S^ from Sj for 
5 = flog(2nD + l)l, 1 < i < D. 

For every i, in a first step, we have S^ = Si and matrices A^ = Dg(si) and 
= Dg~ 1 (si) in ]K rxr . The computation of these matrices does not change the over- 
all complexity. Assume that A^B^ = I mod (t - l) 2m K[[t - 1]]. In the (m + l)th 
step, we evaluate the matrix ^( m + 1 ) = DF(S^ m) ) mod (t - l) 2m+1 K[[t - 1]]. The cost 
of this step is 0(nLi2 m mlog(m)). When m = 0, we take #W = 25W - B^A^B^. 
For m > 1, to define an approximation i?( m+1 ) to the inverse of A^ m+1 \ noticing that 
A ( m +i) = DF{sl m) ) = DFiS^'^) = A^ mod (t - l) 2m_1 K[[t - 1]], we compute 
B' = flW-BW^M-AW^W, We have that A^+^B' = I mod (t - l) 2m K[[t - 
1]]. Then, £( m+1 ) = 2B' - B'A^+^B' is obtained, which satisfies A^ 1 ^^ 1 ) = I 
mod (i — l) 2m+1 ]K[[t — 1]]. The complexity of this step is 0(?i 3 2 m m log(m)). Finally we 
evaluate S^ m+1) = s[ m) - £( m+1 ) F (S^ ') within 0(Li2 m m log(m)) operations. 

The total number of operations over K to do this step for every 1 < i < D is 0(n(nLi + 
n 3 )D 2 log(D) log log(D)). 

Third step: This step consists in the computation of P(0, U, a) = ^2h=oPh{^i Oi)U h and 
$1(0, U, a) = Eh=o fe(0, a)U h for a generic a = (a u . . . , a n ) € Q n . 

By expanding P(t,C7,y) = ELo^ 1 ^ g k [[* " into P owers of ^ (ift " 

ai), . . . , (y n — a n ), we have that the coefficients corresponding to U h and U h (yk — afc) 
(1 < k < n, < h < D) are pj l (t,a)/q(t) and ^-{t,a) / q{t) respectively. As the degrees 
of the polynomials involved in these fractions are bounded by nD (see Lemma [7]), they 
are uniquely determined by their power series expansions modulo (t — l) 2nD+1 ¥L[[t — 1]] 
(see [SI Corollary 5.21]). 

The algorithm proceeds as follows: first, it computes the coefficients of U h and U h (yk — 
a k ) (1 < k < n, < h < D) in P(t, U,y) = Y\Zi( U ~ K^hV)) e K[t][U,y] following gH 
Algorithm 10.3] in 0(n 2 D 2 log 3 (-D) log log 2 (-D)) operations over K. From these coefficients, 
Ph(t,a) and 2j^(t,a) (1 < k < n, < h < D), and q(t) are obtained within complex- 
ity 0(n 2 £> 2 log 2 (.D)loglog(Z))) over K by using [HJ Corollary 5.24 and Algorithm 11.4] 
and converting all rational fractions to a common denominator. Finally, the algorithm 
substitutes t = in these polynomials to obtain P(0, U, a) and ^~(0, U, a) for 1 < k < n. 

Fourth step: The algorithm computes Q(U,a) = gcd(P(0, U, a), §^(0, U, a)) within 
complexity 0(D log 2 (D) log log(D)) and makes the required exact divisions by Q(U,a) 
leading to the geometric resolution. This last step does not change the overall order of 
complexity, which is 0(n 2 D 2 log(D) log log(D) (L + n 2 d + log 2 (£>) log log(D))) . □ 

The algorithm underlying the proof of Proposition[l4]can be adapted straightforwardly 
to handle the cases s = 1 and s = n within the same complexity bounds. 
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4.3 Main algorithm 

The main algorithm of this section is the following: 
Algorithm 15 

Input: Polynomials f\, . . . , f m £ K[x±, . . . , x n ] satisfying Hypothesis encoded by an sip 
of length L, and positive integers d\, . . . , d m such that deg fi < di for 1 < % < m. 

Output: A finite set Ai C A n intersecting the closure of each connected component of the 
realization of every feasible sign condition over fi, . . . , f m encoded by a list Q of geometric 
resolutions of ^-dimensional varieties. 

Procedure: 

1. Make a random linear change of variables with coefficients in Q. 

2. Take a point p = (pi, . . . ,p n ) G Q n at random. 

3. Starting with Q : = 0, for k = 1, — 1 and for every S C {l,...,m} with 
1 < #S <n- k + 1: 

(a) Obtain an sip encoding the polynomials which define the variety Wk t s associated 
with the polynomials . . . ,Pk-i,%k, ■ ■ ■ , x n ), . . . , f m (pi, ■ ■ ■ 

and the projection to the kth coordinate x^. 

(b) Compute a geometric resolution {q( k,s \U), vjf (U), . . . ,v^' S \u)} C K.[U] of 
the variety ir x ( s H {t = 0}) C A n ~ fc+1 by means of a deformation from a type 
1 initial system and add the geometric resolution 

to the list Q . 

4- Add to the list Q the geometric resolutions {fi{pi, ■ ■ ■ ,p n -i, U),p±, . . . ,p n -i, U}, for 
1 < i < m, and {U,p±, . . . ,p n }- 

Theorem 16 Alaorithm\15\ is a probabilistic procedure that, from a family of polynomials 
fx, ■ ■ ■ i fm £ K[xi, . . . , x n ] satisfying Hypothesis obtains a finite set M intersecting 
the closure of each connected component of the realization of every sign condition over 
fi , . . . , f m . If the input polynomials have degrees bounded by d > 2 and are encoded by an 
sip of length L, the algorithm performs 0((^^^ m ' n ^ (™) (™Zi) )d 2n n 4 log(d) (log(n) + 
log log(d)) (L + n 2 + n log 2 ((i)(log(n) + loglog(d)))) operations in K. 

Proof: Assuming that the random linear change of variables made in the first step of the 
algorithm is generic in the sense of Proposition [H by Proposition [3l it suffices to show that 

M U (U U Z(C,k))cM 
k=i cec(fc,p) 
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where C{k,p) denotes the set of all the connected components of the M n -subsets Tn {xi = 
pi, . . . ,Xk-i = Pk-i} with r a connected component of a feasible sign condition over 
fl i ■ ■ ■ j fm ■ 

Note that for a generic point p = (px, • • • ,Pn) £ ^ n > f° r every 2 < k < n, Hypothesis 
[9] also holds for the polynomials fi(pi, ■ ■ ■ ,Pk-i, x k, ■ ■ ■ , %n), I < i < n. Thus, by taking p 
at random, we may assume that the hypothesis is met at each step of the recursion. 

Then, for every 1 < k < n — 1, by Proposition [11] we have that 

(J Z(C,k)c (J {pi,...,Pfc- 1 }x7r x (y fc|S n{t = 0}), 

CeC(k,p) Sc{l,.,m} 
l<#5<n-fc+l 

and Step 3 of the algorithm computes geometric resolutions for the sets in the right-hand 
side union; therefore, UceC(fc p) %(C, k) C M.. Finally, note that UceC(n p) %(C, n) C 
yJiLi{fi(pii ■ ■ ■ >Pn-l,%n) = 0}, which along with the point p, is added to the set M. in 
Step 4 of the algorithm. This proves the correctness of Algorithm [T5l 

For every 1 < k < n — 1 and each S C {1, . . . , m} of cardinality at most n — k + 1, 
the sip encoding the polynomials which define the variety Wk : s computed at Step 2 of the 
algorithm can be taken of length 0(nL+n s ). Moreover, the number of points in V^gR-ft = 
0} is bounded by (™Zi) d n ~ k+1 . Therefore, the result follows using the complexity estimate 
in Proposition 1141 □ 

Now we will show how to get the entire list of feasible sign conditions over the poly- 
nomials fi, . . . , f m satisfying Hypothesis El The algorithm relies on the following: 

Proposition 17 Let fi,---,f m £ M.[x±, . . . ,x n ] be polynomials satisfying Hypothesis^ 
and let M be a finite set such that M n C ^ for every connected component C of the 
realization of each feasible sign condition over fi, ■ ■ ■ , f m - Then, the set of all feasible sign 
conditions over /i,...,/ m is \J a ^c(M)^ >r7 w ^ ere £(.M) is the set of all sign conditions 
satisfied by the elements of M and P a denotes the subset of {<,=, >} m consisting of the 
elements that can be obtained from a by replacing some of its "= " coordinates with "< " 
or ">". 

Proof: Let a be a feasible sign condition and C a connected component of {x S M. n \ 
/i(x)a"i0, . . . , fm(x)& m 0}- Consider a point z G M. n C and let a G £(-M) be the sign 
condition over fi, ■ ■ ■ f m at z. By continuity, it follows that a G P a . 

Now, let a G C(M) and z G M such that fi(z)ai0 for 1 < i < m. Without loss of 
generality, assume a = (=, ...,=,>,...,>) with t "=" and m — t ">". If t = 0, P a = {a}. 
Suppose now t > 0, and let a G P a . We may assume a = (=, ...,=,>,...,>) with / "=", 
where < I < t. Since the vectors V fi(z), . . . ,V ft(z) are linearly independent, there 
exists v G M. n such that (Vfi(z),v) = for 1 < i < I and (Vfi(z),v) > for I + 1 < i < t. 
Consider a C°° curve 7 : [—1, 1] — ► {fi = ■ ■ ■ = fi = 0} such that 7(0) = z and 7'(0) = v. 
For I + 1 < i < t, fi o 7(0) = and (fi o 7)'(0) = (Vfi(z),v) > 0; therefore, for a 
sufficiently small u > 0, /;o 7(7/) > holds. In addition, for 1 < i < I, fi o 7(n) = for 
every u G [—1,1]. Finally, for t + 1 < i < m, as fi o 7(0) > 0, we have fi o *y(u) > for a 
sufficiently small u. We conclude that a is feasible. □ 
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Given a geometric resolution {q(U),Vi(U), . . . , v n (U)} C K[U] consisting of polynomi- 
als of degree bounded by 5, if fx, . . . , f m are encoded by an sip of length L, it is possible 
to obtain the signs they have at the points represented by the geometric resolution within 
complexity 0(L<5 log(<5) log log(<5) + m5 w ') (here u> < 2.376 is a positive real number such 
that for any field k it is possible to invert matrices in k rxr with 0(r t ") operations, see 
[T7]). First, for 1 < i < m, compute fi{v\(U), . . . , v n {U)) mod q(U) within complexity 
0(L5log(5) loglog(<5)) ([1H Chapter 8]) and then, evaluate the signs of these polynomials 
at the roots of q by using the procedure described in \15\ Section 3] within complexity 
0(m<5 w ). 

By applying the previous procedure to the output of Algorithm [15] and using Propo- 
sition [T71 we compute the list of all feasible sign conditions over fi, ■ ■ ■ , f m - 

Theorem 18 There is a probabilistic algorithm that, given polynomials fi, ■ ■ ■ , f m G 
K[xi, . . . , x n ] of degrees bounded by d > 2 satisfying Hypothesis^ and encoded by an sip of 
length L, computes the list of all feasible sign conditions over these polynomials within com- 
plexityO(Y^=i mM (7)((L+n 2 d)(^: i 1 ) 2 d 2 «n 4 log(d)(log(n)+loglog(d)) + m^"Q: i y)). 
□ 

Our algorithms and complexity results can be refined if we are interested in a particular 
sign condition a over /i, . . . , f m : 

Remark 19 Let a G {<,=, >} m and E a = {i | <Tj = " = "}. Due to Proposition [771 in 
the third step of Algorithm [J_5| it suffices to consider those sets S C {1, . . . m} such that 
E a C S. Then, if #E a = I, in the complexities of Theorems{T^ and \18l the sum can be 
taken over I < s < min{m, n} and the combinatorial factor ( m ) can be replaced by ( m ~, ) . 



5 Closed sign conditions over arbitrary polynomials 

In the case of arbitrary polynomials, our approach is similar to the one in [7]. We will 
consider the same kind of deformations as in the previous section but we will use different 
initial systems whose particular properties enable us to recover extremal points from the 
solutions of the deformed systems. 

5.1 Initial systems for deformations 

Let d be an even positive integer and T the Tchebychev polynomial of degree d. 

Definition 20 For a given s G N with 1 < s < n and r := s + n — 1, a type 2 initial 
system is a polynomial system of the form: 

9i{x) = Ti (n + A i{n+1) + Ei<fc<n AkT(x k )^j forl<i<s, 

gi(x,fi)= Hj - 1 — 0) fors + l<i<r, 
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where Tj G {+> — } / or 1 < i < s, and ^4 G Q sx ( n+1 ) ffog Cauchy matrix defined as 
A-ik = a .j)_fc / or l<i<s, l<fc<n + l, with < ai < • • • < a s integers such that 
a s + n + 1 is a prime number. 

For s = 1, a type 2 initial system consists of a polynomial g\{x) as above and its 
partial derivatives d9 g^ , • • • , ^f^~ ■ Finally, for s = n, a type 2 initial system consists of 
n polynomials g\,. . . ,g n constructed as above from the Cauchy matrix A G Q nx ( n + 1 ). 



Note that if n = " + ", then gi{x) > for every x G 
gi{x) < for every x G K n . Moreover, for s + 1 < i < r, 



and, if T{ 



gi(x,fi) = T'(xi- a+1 )( 



T 3 A j(is+l)l l j 
l<j<s 



then 



(11) 



The Bezout number of a type 2 initial system is D 



\)d s {d- l) n - 



< 



Lemma 21 Property (H) (see Section \3.3.1\) holds for any type 2 initial system. 



Proof: Assume 1 < s < n. For s + 1 < i < r, let g~i(x, fii, . . . , /U s _i) = gi(x, fi\, . . . , fi s ~i, !)■ 
Let S C {2, . . . , n} be a set with n—s elements, let e : B — > {—1, 1} and suppose e(k) = 
1 for a elements in B. Let Ss,e be the set of solutions (x, /2) = (xi, . . . , x n , /2i, . . . , fl s -i) 
of the system 

gi (x) = ■■■ = g s (x) = g s+1 (x,/j) = g r (x,n) = 



which also satisfy 



T'(xk) = and T{xk) = e(k) for every k G B. 



(12) 



(13) 



We are going to compute the number of elements of S B, e ■ Without loss of generality, 
suppose B = {s + 1,... ,n}. As gcd(T',T + 1) = T d/2 and gcd(T',T - 1) = T'/T d/2 
(where T d / 2 is the Tchebychev polynomial of degree d/2), the number of (n — s)-uples 
(x s+ \, . . . ,x n ) satisfying (fT3|) is (d/2) n ~ s ~ a (d/2 — l) a . Each of these elements can be 
extended to solutions (x,p,) of (|12p in d s different ways: Let A' G <Q SXS be the Cauchy 
matrix formed by the first s columns in A. The conditions gi(x) = ■ ■ ■ = g s {x) = and 
PD, imply 

( T ^i)\ ( n + A 1{n+1) + A lk e(k)\ 

s+l<k<n 



A' 



\ T(x s ) J 



n + A s(re+1) + ^ A sk e(k) 

s+l<k<n 



(14) 



Using Cramer's rule for the previous system, we deduce that, for each solution (xi, . . . , x s ), 
T{x~k Q ) 7^ il (1 < ^ s ) because all these numbers are rational with denominator a 
multiple of the prime number a s + n + 1 and numerator relatively prime to it. As the 
equation T{x) = a only has multiple solutions for a = ±1, for every 1 < fco < s, there 
are exactly d values X£ for which T(xk ) satisfy (|14p . It can be easily shown that each 
solution (xi, . . . , x n ) of (fl~3l) and (fT4|) can be extended uniquely to a solution (x, p) of (fT2l) , 
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As for (x,a) G SB, e , T(xk) = e(k) = ±1 for every k G B, the sets Sb,b are mutually 
disjoint. Then, if (iW.jjW), (x( 2 ),/l( 2 )) are two different solutions of (|12p . x^ 1 ) 7^ x^ 2 ). By 
taking into account every B C {2, ... , rt}, every a (0 < a < n — s) and every function 
e : -B — > {—1,1} whose value is 1 at exactly a elements in B, we find 



E 

0<a<n— s 



n — s 
a 



dsn-s-a ,d 
2J V2 ~ 



(d-l) n ~ S d S 



solutions to (fT2j) . 

Consider now the Jacobian matrix of this system evaluated at each of these solutions 
and suppose, without loss of generality, that the solution (x, p) considered corresponds to 
B = {s + 1, . . . , n}. Then, this matrix is of the form 



* { 










s-1 { 


{: 





C2 


n — s { 




c 3 


* 



with C±, C2 and C3 invertible matrices. 

For s = 1 and s = n, the proof is similar. 



s-1 



□ 



5.2 Geometric properties 

Let /1, • • • , / m G R[xi, . . . , x n ] and let S := {ii, . . . , i s } C {1, . . . , m} with 1 < s < n. 
As explained in Section 13.21 for every point z G M n with maximum or minimum first 
coordinate over the set {x G ffi n | /^(x) = 0, . . . , fi s (x) = 0}, there exists a nonzero vector 
a = (/ii, . . . , fi s ) such that (z, fx) is a solution to the system 



0/<, 



(15) 



I ^i<i<« 9x2 ( z ) — " " — Si<j< s &^ ( z ) — 0- 

Now, a homotopic deformation of this system by means of a type 2 initial system is as 
follows: for every 1 < k < s, Fk(t, x) = (1 — t)fi k + tgk(x), and, for every s + 1 < k < r, 

dfi, 



1 



dx 



fc-s+l 



(x) + * Z 



<9x 



fc-s+l 



E 



9^ 



<9x fc 



A--.S+1 



■(t,x). 



Thus, for any to G R and every xq G R n at which the first coordinate function attains a 
local maximum or minimum over the set {x G R n \ Fi(to, x) = ■ ■ ■ = F s (to, x) = 0}, by the 
implicit function theorem, there is a nonzero vector uq G R" such that Fi(to,xo) = ■ ■ ■ = 
F s (t ,x ) = F s+1 (t ,x ,fi ) = ■■■ = F s (t Q ,x Q ,u ) = 0. 

In the sequel we will consider deformations by means of specific type 2 initial systems. 
Let d G N be an even positive integer with d > deg/j for every 1 < i < m. Let 
qx < ■ ■ ■ < q m be the first m prime numbers greater than n. For 1 < i < m, let 

1 



Si 0*0 



1 \ - 

n+ - + 



Kfc<n 



n 



1 + fe 



T(x k ) and ^ (x) = -gMx) 
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Note that for each S = . . . , i s } C {1, . . . , m} with 1 < s < n and every list r±, . . . , t s of 
+ and — signs, the polynomials gj 1 , . . . , gj s form a type 2 initial system with a,j = qi j —n— 1 
for 1 < j < s (see Definition I20p . In addition, for 1 < i < m, we denote 

F+(t,x) = (l-t)fi(x)+tgf(x) and Fr(t, x) = (1 - t)fi(x) + t g r(x). 

Lemma 22 Let S = {ii, . . . , i s } C {1, . . . , m} with s > n and n, . . . , r s a list of + and 
— signs. Then, the set {t G C | 3x G C n with i^ Tl (t,a;) = • • • = F^ B (t,x) = 0} is finite 
(possibly empty). 

Proof: Denote by F?} , . . . , FJ S , gj 1 , . . . , gT s the polynomials obtained by homogenizing 
F^ , . . . , F[ s s and gj^ , . . . , gj* with a new variable x$. 

Let Z C A 1 x F n be the set of common zeros of F£ , . . . , £[* and vr t : A 1 x P n -> A 1 the 
projection. In order to prove the statement it suffices to show that nt(Z) is a finite set. 
Since 7r t is a closed map, this can be proved by showing that 1 ^ 7it(Z), or equivalently, 
that the system gV^x) = ■ ■ ■ = g1 s {x) = has no solution in P n . 

First, note that, if (1 : x\ : ■ ■ ■ : x n ) is a solution to this system, then (T(xi), . . . , T(x n )) 
is a solution to the linear system B.y 1 = — (n+ . . . ,n + — — — )*, where B S Q( n + 1 ) xn is 

the Cauchy matrix of coefficients of gj* , ... , g^l ■ But this linear system has no solutions, 
since its augmented matrix has a nonzero determinant (as in the proof of Lemma [2TI it can 
be shown that this determinant is a rational number whose numerator is not a multiple 

Finally, we have that, for 1 < j < s, g> (0, sci, . . . , x n ) = ±2 d_1 J2i<k<n a -n-i+k x i- 

J — — Hlj T 

Considering the equations for 1 < j < n, we deduce that (xf, . . . , is in the kernel of the 
Cauchy matrix ( — _ n 1 _ 1+fc )i<j,fc<n an d therefore, it is the zero vector. Thus, the system 

gj^(x) = ■ ■ ■ = gj°(x) = has no solutions in {xq = 0}. □ 

Notation 23 For S = . . . , i s } C {1, . . . , m} with 1 < s < n, and r = (n, . . . , r s ) E 

{+, —} s , we denote Vs, T C A 1 xA n xP s ~ 1 i/ie variety defined by the polynomials constructed 
as in (0) by taking hi, ■ ■ ■ , h r as the polynomials in system U5\) and gi,...,g r the type 2 
initial system given by g[^ ,...,gT\ We consider the decomposition Vs,t = Vg°J U V^ 1 ^ U Vs, T 
as in (0). 

The following proposition will enable us to adapt Algorithm [15] in order to solve the 
problem in this general setting. 

Proposition 24 Let a G {<,=,>} m , E a = {% \ Oi = " = "}, U a = {i \ <n = " > "} 
and L a = {i \ o~i = " < "}. For S = {ii,...,i s } C {1, ...,m} with 1 < s < n, 
let T s = {t e | Tj = " + " z/ij G [7 CT and Tj = " - " i/i^ G L a }. Then, 

after a generic linear change of variables, for each connected component C of {x G 
W 1 | fi(x)ai0, ... , f m (x)cr m 0}, we have 

z (°) c U U *x(v s , T n{t = o}). 

SC{l,...,m} tSTs 

l<#S<n 
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Proof: Without loss of generality, we may assume that E a = {1, . . . , I}, U a = {l+l, . ■ ■ , k}, 
and L a = . . . , m} for some /, k with < I < k < m; that is, we consider a connected 

component C of V = {/i = 0, . . . , ft = 0, fi +1 > 0, . . . , f k > 0, / fc+ i < 0, . . . , f m < 0}. By 
Proposition [H after a generic linear change of variables, Z(C) is finite. Moreover, since V 
is a closed set, Z(C) C C. Let z G Zi n f(C) and < e < 1 such that: 

• B(z,e) n P C C and n Z(C) = {z}, 

• for every S = . . . , ig} C {1, . . . , m} with s > n and every (n, . . . , Tg) G {+, — } s , 
e < |to| for every t in {t G C \ {0} | 3x G C n with F£ (t, ac) = • • • = F[° (i, x) = 0}, 

• for every S C {1, . . . , m} with 1 < #S < n and every r G 7s 5 £ < |*o| f° r every 
to G C such that has an irreducible component contained in {t = to}- 

For t G R, let R t be the set 

{x G B(z,e) | F x + (t,x) > 0,...,F+(i,x) > 0,Ff (t,z) < 0, . . . , F~(t, x) < 0, 
F+ 1 (t,x)>0,...,F+(t,x)>0,F- +1 (t,x)<0,...,F-(t,x)<0}. 

We have that Rq = C n e) and, for every t g [0,1], z 6 Rt- 

Let > be the distance between the compact sets dB(z,e) n {xi < z\] and i?o. 

We claim that for some t\, < t\ < e, the connected component C 1 of containing z 
is included in {x G B(z,e) \ cL(x,Rq) < v/2}. Suppose this is not the case. Let (^) n eN De 
a decreasing sequence of positive numbers converging to and with t'-y < e, and for every 
n G N, let C' n be the connected component of R t > containing z. Since C' n intersects {x G 
B(z,e) | d(x,Ro) > v/2}, there is a point r n G such that d{r n ,Ro) = v/2. Consider 
the sequence (r n ) nS N C B(z,e), and a subsequence (r rej .)j 6 N convergent to an element 
r G B(z,e); then, d(r,Ro) = v/2. Now, for 1 < z < k, we have that F^(t nj ,r nj ) > 
for every j G N, and so, i^ + (0, r) > 0, and for every 1 < i < I and k + 1 < i < m, 
Fi~ {tn 3 •) ) < 0, which implies that F~(0, r) < 0. Therefore r G -Ro> contradicting the 
fact that d(r, R Q ) = v/2 > 0. 

Let w G C be a point at which the function xi attains its minimum over C . Since 
z G C, we have it?i < z±. If w 6 8B(z,e), then w G dB(z,e) n {xi < zi}, and so, 
d(w,Ro) > v, contradicting the fact that d(w,Ro) < v/2. Therefore, w G B(z,e). 

As each of the polynomials F± , . . . , F^ , F± , . . . , F { ~ , , . . . , F^~ , i 7 ^.]^ , • • • , F^ that 
does not vanish at (ti,w) takes a constant sign in a neighborhood of this point, we conclude 
that, if F?},... } F[* are all the polynomials vanishing at (t\,w), then the function x\ 
attains a local minimum over the set {x G R n | F^(t±,x) = 0, . . . , F^ s (ti, x) = 0} at 
it?. Let 5o = . . . ,i s }, which is not empty. For 1 < i < m, F^(tx,w) and F~(ti,w) 
cannot be both zero; then, i\, . . . ,i s are all distinct. Because of the way we chose s, we 
also have that s < n. Now, if to = (ti,...,t s ), we have that (ti,w) G vr ii2; (Vs 0>T0 ), but 
taking into account that < t\ < s, it follows that (ti,w) G Kt,x(Vs ,T )- Therefore, 
(h,w) G Usc{i,.,m} U T eT q *t,x( v S,T) and < \(h,w) - (0,z)\ < y/2e. 

l<#S<n S 

Since the previous construction can be done for every e > sufficiently small and the 
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sets TTtx{Vsr) are closed, we conclude that (0, z) £ (Jsc{i,...,m} UtpT zC^St)- Therefore, 

l<#S<n S ' ' 

z£ U U T*(Vs,rn{t = o}). 

Sc{i,...,m} reT s 

l<#S<n 

□ 

5.3 Symbolic deformation algorithm 

In the sequel, £1 will denote a positive real number such that for any ring R, addition, 
multiplication and the computation of determinant and adjoint of matrices in R kxk can 
be performed within 0(k : ) operations in R. We may assume Q < 4 (see [11]) and, in 
order to simplify complexity estimations, we will also assume that Q > 3. 

Proposition 25 There is a probabilistic algorithm that, taking as input polynomials hi,... ,h r 
in K[xi, . . . ,x n ,fii, . . . ,fjb s ] as in encoded by an sip of length L, obtains a geometric 
resolution of a finite set containing ir x (Vr\ {t = 0}) for a deformation defined from a type 
2 initial system within complexity 0(n 3 (L + dn + n n ~ 1 )D 2 log 2 (D) log log 2 (£))) , where d 
is an even integer such that d > deg x (hi) for every 1 < i < r. 

Proof: The structure of the algorithm is similar to that of the algorithm underlying the 
proof of Proposition [T5J 

First step: Take a = (ai, . . . ,a n ) £ Q n at random and compute a geometric resolution 
associated to the linear form £ a (x) = a±x± + ■ ■ ■ + a n x n of the variety defined in A r by 
the (dehomogenized) type 2 initial system. 

As shown in the proof of Lemma [2T1 this variety can be partitioned into subsets Sb^- 
So, we first compute a geometric resolution associated with £ a (x) for each Ss,e'- after solv- 
ing a linear system of the type (fl4l) . the x-coordinates of points in Ss,e turn to be defined 
by a square polynomial system in separated variables; then, the required computation can 
be achieved as in [271 Section 5.2.1] within complexity 0(D% e log 2 (D B)e ) loglog(Z) S) e)), 
where Ds, e is the cardinality of Ss,e- 

Finally, a geometric resolution of the whole variety is obtained following the split- 
ting strategy given in [441 Algorithm 10.3], and noticing that, if {q, qo, w±, . . . , w r } and 
{q, q~o, wi, . . . , w r } are geometric resolutions of disjoint sets with q and q coprime poly- 
nomials, then {qq, q$q + %q, w\q + Wiq, ■ ■ ■ , w r q + w r q} is a geometric resolution of their 
union. This can be done within 0{nD log 2 (Z?) log log(D)) operations in Q. 

The whole complexity of this step is 0(nD 2 log 2 (D) log log(D)). 

Second step: Compute P(t, U, y) mod ((* - l) 2 " D + 1 + ( yi - a± , . . . , y n - a n ) 2 )K[[t - 

mu,y}. 

First, from the geometric resolution computed in the previous step, obtain a geometric 
resolution associated with £(x, y) = y±xi + • • • + y n x n of the variety defined by the initial 
system over modulo the ideal (yi — a>i, . . . ,y n — a n ) 2 , applying [20l Algorithm 1] 

within complexity 0((dn 2 +n^)D 2 log 2 (D) log log 2 (D)). Then, consider the variety defined 
by Fi, . . . , F r over K(t, y) (see Lemma [6]). Since iq, . . . , F r can be encoded by an sip of 
length L + 0((d + s)n), a geometric resolution of this variety associated with the linear 
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form £(x, y) modulo the ideal (t— l) 2nD+1 + (yi — ct\, . . . , y n ~ a n) 2 can be obtained from the 
previously computed geometric resolution by applying \20\ Algorithm 1] within complexity 
0(n 3 (L + dn + n n ~ l )D 2 \og 2 {D) log log 2 (L>)). 

Third step: From the approximation to P(t, U,y) obtain the required geometric resolu- 
tion, by performing the same computations as in the third and forth steps of the algorithm 
underlying the proof of Proposition [141 which does not modify the overall complexity. □ 

The algorithm underlying the proof of Proposition[25lcan be adapted straightforwardly 
to handle the cases s = 1 and s = n within the same complexity bounds. 

5.4 Main algorithm 

Here we prove the main result of this section. 

Theorem 26 Given polynomials /i, . . . , f m £ M.[x\, . . . , x n ] with degrees bounded by an 
even integer d and encoded by an sip of length L, for generic choices of the param- 
eters required at intermediate steps, the algorithm obtained from Algorithm [75| taking 
p = (0, . . . , 0) and replacing step 3. (b) with 

(b ! ) For every r £ {+, — , compute a geometric resolution {q^ k,s,T \U),v^' S ' T \u), ... , 

v£' S ' t) (U)} C K[U] of a finite set containing ir x (V k)S ,r n {t = 0})) C A n ~ fc+1 by 
means of a deformation from a type 2 initial system, and add the geometric resolution 

{q^\U), 0, . . . , 0, v[ k ^\u), vt S ' T \U)} 

to the list Q. 

computes a finite set Ai C A n intersecting each connected component of the realization 
of every feasible closed sign condition over fi, . . . , f m . The complexity of the algorithm is 
0(n*(L + d+ n 2 ) log 2 (d)(log(n) + loglog(d)) V«( 2 *(™) (^) 2 )). ^ 

Proof: As in the proof of Theorem [T6l by Proposition [31 it suffices to show that 
Ui<fc<n-i UceC(fcp) Z(C, k) C A4, which is a consequence of Proposition [2^1 

Taking into account the linear change of variables performed at the first step of the 
algorithm, for 1 < k < n — 1 and every S C {1, . . . ,m} with 1 < #S < n — k + 1, 
we can obtain an sip of length 0(nL + n 3 ) encoding the polynomials involved in system 
(|15|) . Then, the algorithm underlying the proof of Proposition 1251 computes a geometric 
resolution of a finite set containing ^(T^; St H {t = 0})) C A n ~ The stated complexity 
is obtained by adding up the complexities of these steps for all (k, S,r). □ 

Proceeding as in Theorem 1181 from a finite set A4 intersecting each feasible sign con- 
dition over the polynomials fi, . . . , f m , we can obtain the list of all closed sign conditions 
over these polynomials by evaluating their signs at the points in M. We deduce the 
following complexity result: 
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Theorem 27 There is a probabilistic algorithm which, given polynomials /i,---,/ m £ 
. . . , cc n ] of degrees bounded by an even integer d and encoded by an sip of length 
L, computes the list of all feasible closed sign conditions over these polynomials within 
0(ES {m ' n} (™)(:: i 1 ) 2 d 2re (2 s n 6 (L+dW)log 2 (d)(log(n)+loglog(d))Vm^- 2 )-(^ 1 1 )" 
operations in K. □ 



6 Some particular cases 
6.1 The bivariate case 

Here we will show that when n = 2, Algorithm 1151 solves our main problem for an arbitrary 
finite family of polynomials. 

Lemma 28 Let f G C[aci,X2] be a nonzero polynomial with no factors in C[xi] \ {0} 
and f = ni<-;<a^ irreducible factorization in C[xi,x 2 \. Let g\,g 2 be polynomials 
satisfying conditions (H), with g\ relatively prime to f . Let F\ = (1 — t)f + tgi and 
F 2 = (1 — t)fa^ + tg2 and V be the variety defined in $8$). If z G C 2 satisfies that either 

there is an index Iq such that Pi (z) = -^r-(z) = or there are indices i\ ^ i 2 such that 
p h (z) = p i2 (z) = 0, then z G ir x (V n {t = 0}) . 

Proof: In order to simplify the notation, we write /' = J^-. Set / = (F\, F 2 ) C K.[t, x\, x 2 ]. 
Then V n {t = 0} = V((I : t°°) + (t)). 

Let h = gcd(/,/') = nip? -1 , = ///» = TliPi and ^ 2 = f'/h = >l, .. ; !>/)■ 
We will show that the equalities (/ : t°°) = (I : t) = (F 1 ,F 2 ,h 2 gi - ^132) hold. Then, 
since each condition in the lemma implies that h\(z) = h2{z) = 0, we deduce that (0, z) £ 
V(I : t°°) and, therefore, that z G 7r z (y n {t = 0}). 

To prove the second equality, first note that (/i2<?i — hig 2 )t = h 2 F\ — h±F 2 , which 
shows the inclusion D. Now, if p(t,x) G (/ : t), we have p(x,t)t = (a.\{t,x)t + ao(x))F\ + 
(/3i(t,x)t + /3q(x))F 2 for polynomials «i, ao, /3o- Substituting t = 0, we obtain aof = 
—0of' and, dividing by h, a.Qh\ = —(3oh 2 . Then, there exists c G C[x] such that ao = ch 2 
and /3o = —ch\ and therefore, p(t,x) = a\{t,x)F\ + Pi(t,x)F 2 + c(h 2 g\ — h\g 2 ). 

The first equality will be proved by showing that (I : t 2 ) C (I : t). Let t) G (I : t 2 ). 
Then, p(x,t)t 2 = (a 2 (x,t)t 2 + a\(x)t + ao(x))Fi + (b 2 (x,t)t 2 + bi(x)t + bo(x))F 2 for certain 
polynomials ai,bj. Setting i = in this identity, it follows that there is a polynomial 
c\ G C[x] such that ao = c\h 2 and 60 = —c\h\. Now, by looking at the terms of degree 1 
in t, we have that a$gi + bog 2 = —a±f — b±f. Then, c\{h 2 g\ — h\g 2 ) = h{—a\h\ — b\h 2 ). 
We deduce that h divides c±: for every i with Si > 1, we have that pi \ h\ and p. t \ h 2 g\ 
(since / and g\ are relatively prime); so, pi \ (h 2 g\ — h\g 2 ) and, therefore | c. If 

c 2 G C[x] satisfies c\ = c 2 h, we have p(x, t)t = (a 2 t + ai)i*i + (62* + &i)-^2 + c 2 (f'gi — fg 2 )- 
Since /'oi - fg 2 = -{g 2 - f')Fi + (gi - f)F 2 G J, we conclude that p £ (I : t). □ 

Proposition 29 Let fi,...,f m be arbitrary bivariate real polynomials and a G {<,=,> 
} m . Then, after a generic linear change of variables, for each connected component C of 
{x G W 1 I /i(x)o-i0, . . . , / m (x)cr m 0} we have that Z{C) C Usc{l,...,m}, 1<#S<2 ^(Vs fl {t = 
0}) where the varieties Vs are defined from type 1 initial systems. 
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Proof: After a generic linear change of variables we may assume that, for each connected 
component C, either nx(C) =1 or Z(C) is a non-empty finite set (see Proposition [1]). 

Assume Z m {(C) / and let z = (^1,^2) G -^mf(C). Since z G dC, there is an index 
io such that fi ^ and fi (z) = 0. If /j has two or more non-associate irreducible 
factors in C[x] vanishing at z, or an irreducible factor vanishing at z and whose derivative 
with respect to x<i also vanishes at z, by Lemma [28l z G ^x(Vn \ n {t = 0}) (note that, 
because of the generic change of variables, fi does not have factors of the form X{ — a). 
Otherwise, there is a unique irreducible factor p of fi vanishing at z which must have all 
real coefficients (since its complex conjugate also divides /j and vanishes at z) such that 

fi(*)*o. 

By the implicit function theorem applied to p at the point z, there is a continuous 
curve (xi,X2(xi)) defined in a neighborhood of z\, and a neighborhood of z such that the 
polynomial p (as well as any power of p and also fi ) has constant signs above, below and 
on the curve in this neighborhood. Since Z\ = inf 7Ti(C), there must be an index i\ 7^ io 
such that fi x (z) = 0. Moreover, we may assume that has a unique irreducible factor q 
vanishing at z that is not an associate to p. In this second case, z is an isolated point of 
W {io,h} = V (fioifh) and then > b y LemmaEl z G ^(V^^} n {t = 0}). □ 

Using this proposition, following the proof of Theorem [TBI we have: 

Theorem 30 Alaorithm\15\is a probabilistic procedure that, given polynomials fi, ■ ■ ■ ,f m £ 
K[xi,X2] of degrees bounded by d > 2 that are encoded by an sip of length L, obtains a 
finite set A4 intersecting the closure of each connected component of the realization of 
every sign condition on the polynomials within complexity 0(m 2 d 4: log(d) log \og(d)(L + 
log 2 (d)loglog(d))). □ 



6.2 A single polynomial 

In this section we will show the procedure described in Section [5] can be adapted to solve 
the problem of computing a point in the closure of each connected component of {/ = 0}, 
{/ > 0} and {/ < 0} for an arbitrary polynomial / G K[x±, . . . , x n ]. 

Let d be an even positive integer such that d > deg(/). Let T be the Tchebychev 
polynomial of degree d. We define g(x) = n + - + Ylk=i q - n -i+k T( x k), where q is the 
smallest prime greater than n; F(t,x) = (1 — t)f(x) + tg(x), and, for 2 < i < n, Fi(t,x) = 

Note that g > over R n and, therefore, F > over the set {/ = 0}. Moreover, this is 
a deformation with a type 2 initial system (see Definition |20|) . As in Section 13.3.11 we let 
V = {F = F2 = ■ ■ ■ = F n = 0} C A 1 x A™, and consider its decomposition 

Lemma 31 After a generic linear change of variables, for each connected component C 
of {/ = 0}, {/ > 0} or {/ < 0}, we have Z{C) C tt x (V n {t = 0}). 

Proof: Consider first a connected component C of {/ = 0} (for a similar approach in this 
case with an alternative deformation, see [S])- Let z G Z ia ^{C). Take e > such that 
B(z,e) meets neither a connected component of {/ = 0} different from C nor the finite 
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set Z(C) \ {z}, and e < \to\ for each £ C such that has a connected component 
contained in {t = to}. 

Let [i G (zi, zi + e) be such that dB(z, e) n C C {xi > //}. Without loss of generality, 
we may assume that / is positive over the compact set dB(z, e) n {xi < fi}. Let eo G (0, e) 
be such that F is positive over [— £o,£o] x (9B(z,e) n {xi < //}). Now, let y G B(z,e) 
with yi < zi (thus, /(y) / 0). Since F(— e , z) < 0, F(e , z) > and F(0, y) / 0, there is 
a point (ti,y) in the union of the line segments joining (— £q,z), (0,y), and (0, y), (eojz) 
respectively, such that F{t\,y) = 0. We have ii ^ and y*i < Zi. Let w G {x G B(z,e) \ 
F(ti,x) = 0} be a point at which the coordinate function xi attains its minimum over 
this compact set. Note that w G" dB(z,e), since w\ < y\ < Z\ < [i and F is positive over 
[— £0)£o] x (9B(z, e) n {xi < /i}). Therefore, u; G B(z,e) and so, G V'. Moreover, as 

< \t\\ < e, we have (ti,w) G V. Since < \(ti,w) — (0, z)| < V2e, and the construction 
can be done for an arbitrary sufficiently small e > 0, it follows that (0, z) G V. 

Assume now that C is a connected component of {/ > 0} and let z G Z- m f{C) (then, 
f(z) = 0). Let C be the connected component of {/ = 0} containing z, and e > such 
that B(z,e) meets neither a connected component of {/ = 0} different from C nor the 
finite set Z(C) \ {z}, and e < \to\ for every to G C such that has an irreducible 
component included in {t = to}. 

Let ji G {z\,z\ + e) be such that dB(z,e) nC C {xi > ^i} and 7 : [0,1] — > R" a 
continuous semialgebraic curve such that 7(0) = z and 7((0, 1]) C C H -B(z, e) n {xi < //}. 
Let Ci be the connected component of C n S(z, e) with 7((0, 1]) C Ci. Take t\ G (— e, 0) 
small enough so that F{t\,^{l)) > 0. Since F(£i,7(0)) < 0, there exists u G (0,1) such 
that F(ti,-y(u)) = 0. Let C be the connected component of {x G B(z,e) \ F(ti,x) = 0} 
containing 7(u). As 7(u) G C fl Ci, C" U Ci is a connected set. Therefore C C Ci, as 
C" C -B(z, e)\C and Ci is a connected component of this set. Now let K = C' U (-B(z, e) n 
{xi > fj}), which is a compact set, since C 7 = C U (dB(z,e) DC 7 ) C C" U (dB(z, e) n C) C 
-fT. If to G K is a point at which the function x\ attains its minimum over K, then 
w B(z, e) fl {xi > /x}. Then, to is a minimum of xi over the set C n -B(z, e) n {xi < fj,}. 
Therefore, (ii, to) G V", but since < |ii| < e, we have (t±,w) G V. As before, we conclude 
that (0, z) G V. □ 

Using this lemma, due to Proposition [3] and proceeding as in the proof of Theorem [26l 
we have: 

Theorem 32 Given a polynomial f G K[xi, . . . , x n ] of degree bounded by an even integer 
d and encoded by an sip of length L, for generic choices of the parameters required at 
intermediate steps, the algorithm obtained from Algorithm] 15\ modifying step 3.(b) in order 
to deal with type 2 initial systems, computes a finite set M. C A n intersecting the closure 
of each connected component of the sets {/ < 0}, {/ = 0} and {/ > 0}. The complexity 
of the algorithm is 0(n 5 (L + nd + n c_1 ) log 2 (ef) ( log(n) + loglog(d)) 2 (i 2n ). □ 
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